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ABSTRACT 

In this paper we explore several structural and dynamical effects on the projected 
velocity dispersion as possible causes of the fundamental plane (FP) tilt of elliptical 
galaxies. Specifically, we determine the size of the systematic trend along the FP in the 
orbital radial anisotropy, in the dark matter (DM) content and distribution relative to 
the bright matter, and in the shape of the light profile that would be needed to produce 
the tilt, under the assumption of a constant stellar mass to light ratio. Spherical, non 
rotating, two-components models are constructed, where the light profiles resemble 
the R 1 ' 4 law. For the investigated models anisotropy cannot play a major role in 
causing the tilt, while a systematic increase in the DM content and/or concentration 
may formally produce it. Also a suitable variation of the shape of the light profile can 
produce the desired effect, and there may be some observational hints supporting this 
possibility. However, fine tuning is always required in order to produce the tilt, while 
preserving the tightness of the galaxies distribution about the FP. 
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1 INTRODUCTION 

Elliptical galaxies do not populate uniformly the three di- 
mensional parameter space having as coordinates the central 
velocity dispersion a a , the effective radius R c , and the mean 
effective surface brightness I c — Lb/2-kRI, where Lb is the 
total galaxy luminosity in the blue band. They rather closely 
cluster around a plane (Dressier et al. 1987, hereafter D87; 
Djorgovski & Davis 1987; Bender, Burstein & Faber 1992, 
hereafter BBF; Djorgovski & Santiago 1993; and references 
therein) thus called the Fundamental Plane (FP). The exis- 
tence of these "scaling relations" is believed to be of great 
importance for several reasons, including the understand- 
ing of the formation and evolution of elliptical galaxies, and 
their use as tracers of bulk motions, and potentially as a cos- 
mological probe when studying the FP relations for clusters 
at higher and higher redshift. 

For their sample of Virgo ellipticals, BBF introduced 
a convenient coordinate system, where the axis are linear 
combinations of the observables log cr 2 , log-R e and log/ e : 



(1) 



fci = (log o\ + log R, 
fc 2 = (log al + 2 log / c - log R a )/V~6 
fc 3 = (log o 2 - log I c - log R c )/ a/3, 



and the FP is seen edge-on when projected on the fci — &3 
plane. The FP for Virgo ellipticals is shown in Fig. 1 and 



follows the relation: 



k 3 = 0.15 fci +0.36, 



(2) 



having assumed a Virgo distance of 20.7 Mpc, and measuring 
(To, -Rc, and J c respectively in kms -1 , Kpc, and Lb©PC -2 
units (cfr. BBF). The two main properties of the FP for 
Virgo ellipticals are the so called tilt, i.e., the systematic 
increase of k$ along the FP described by equation (2) , and its 
tightness, i.e., the nearly constant and very small dispersion 
of &3 at every location on the FP, with a(k 3 ) ~ 0.05. 

Using the virial theorem, the k-s can be related to the 
total galaxy mass M = C2-R c o" 2 by: 
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(3) 



(4) 



If the virial coefficient C2 is constant for all the galaxies, 
the observed FP tilt as described by equation (2) implies a 
systematic trend in the mass to light ratio with galaxy lu- 
minosity: A//Z/B oc Lb ' 2 (e.g., D87). Meanwhile, the small 
and constant thickness of the distribution about the FP cor- 
responds to a very small (<12 per cent) dispersion of M/Lb 
for any given luminosity. The "smallness" of the 0.2 expo- 
nent may give the impression of the FP tilt being just a 
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minor, not especially important effect. Yet, this is a mislead- 
ing impression. As galaxies in Fig. 1 span a factor ~ 200 in 
luminosity, the tilt corresponds to a factor ~ 3 increase of 
M/Lb along the FP, from faint to bright galaxies. 

On the other hand, the virial coefficient C2 depends on 
the mass, light and velocity dispersion distributions within 
the galaxy, and is constant only insofar all such distribu- 
tions are homologous. It can be obtained by solving the 
Jeans equations, for any reasonable assumption concerning 
the distribution of the bright and DM components, and the 
orbital anisotropy, with the additional constraint that the 
star density distribution will reproduce the observed galaxy 
surface brightness profile. 

Considering only the stellar component (i.e., the frac- 
tion of the total mass whose density follows the light distri- 
bution), its mass M» can be expressed as: 

M t =c* 2 R c a 2 , (5) 

where c\ = C2 in the case of a galaxy that is devoid of dark 
matter (M* = M), while and in general: 

M. M 

— = — • (6) 

Substituting (6) in (3) and (4) we see that the origin of 
the FP tilt can be sought in two orthogonal directions: 
either the tilt may arise from a stellar population effect 
(T* = M»/Lb oc Lb ' 2 while c 2 =const), or from a struc- 
tural/dynamical effect (C2 oc Lb while Y„=const). In 
the former case, the tilt would result from a trend in some 
combination of typical stellar metallicity, age, and initial 
mass function (IMF). A priori this appears to be a quite vi- 
able options: after all, a systematic trend in colours and 
line strengths is known to exists with galaxy luminosity 
(hence with fci), which is usually ascribed to a trend in the 
mean metallicity with the depth of the galactic potential 
well. However, the metallicity effect has been estimated to 
be marginal (D87; Djorgovski & Santiago 1993), and indeed 
existing population synthesys models appear to fail to repro- 
duce but a small fraction of the tilt, unless special conditions 
are verified (Renzini 1995). On the other hand, a drastic 
variation of the IMF along the FP is required to produce 
the tilt, with M/brown dwarfs turning from being a minor 
constituent to dominate the baryonic mass of ellipticals, yet 
with a very small dispersion in the IMF at any location on 
the FP (Renzini & Ciotti 1993, hereafter RC). Searching 
for the origin of the FP tilt in this direction will be further 
pursued in a separate paper (Maraston, Renzini & Ritossa 
1996). 

In this paper we concentrate instead on the second op- 
tion, assuming a constant stellar mass to light ratio Y„, and 
exploring under which conditions structural/dynamical ef- 
fects may cause the tilt in ks via a systematic decrease of 
c|. Although a contribution to the tilt may derive also from 
a trend in the rotational support (decreasing from faint to 
bright ellipticals; e.g., Davies et al. 1983), we concentrate 
here on the effects of systematic trends along the FP in 
1) the degree of radial anisotropy of the velocity dispersion 
tensor, 2) the DM fraction and/or distribution within the 
galaxies, and 3) the density profile of the bright component. 
Concerning the dark matter, a preliminary exploration lead 
to conclude that the central regions of ellipticals should turn 
from baryon dominated to DM dominated with increasing 



Figure 1. The distribution of Virgo (closed boxes) and Coma 
(crosses) ellipticals in the (ki, &2, £3) space, from BBF. The upper 
panel shows the FP edge-on; in the lower panel the FP is seen 
nearly face-on. Open circles represent our reference models Gi 
(see Section 4.5). 

Lb, again with fine tuning required to account for the tilt 
and yet preserve the observed small and constant thickness 
of the FP (RC). A break in the structural homology as a 
possible origin of the tilt has been suggested by Djorgovski 
(1995) and Hjorth & Madsen (1995). 

In Section 2 we briefly describe the dynamical models 
that we have used for exploring points 1) and 2), and derive 
an analytical approximation for the virial coefficient C2. In 
Section 3 we investigate anisotropy as a possible cause of 
the FP tilt, while in Section 4 we explore the effects of the 
amount and distribution of DM, in every case having as- 
sumed a fixed shape of the light profile for all the galaxies. 
Projected velocity dispersion profiles are computed for mod- 
els that succeed in producing the k$ tilt. In conjunction with 
available observations these profiles are then used to reject 
some classes of models and to suggest a future observational 
test for those that still survive. Having completed the anal- 
ysis of the dynamical options for the origin of the tilt, in 
Section 5 we pass to investigate the morphological option. 
This is accomplished by assuming isotropic models without 
DM, whose surface brightness distribution is described by 
R 1 '™ profiles, thus ascribing to a systematic variation of m 
the origin of the tilt. Finally, in Section 6 we discuss and 
summarize the results. 

2 THE MODEL GALAXIES 

We make use of four classes of two-components, spherical 
galaxy models, totally velocity-dispersion supported. We in- 
dicate with r the spatial radial coordinate and with R the 
projected one. 



The Fundamental Plane of Galaxies 3 



2.1 Bright and dark matter distributions 

As well known, the empirical R 1 ^ 4 law (de Vaucouleurs 1948) 
suitably fits the observed surface brightness profiles of many 
elliptical galaxies, but its deprojection cannot be expressed 
in terms of elementary functions. For this reason we describe 
the stellar component by two density distributions which 
give a good approximation to the R 1 ^ 4 law when projected, 
and at the same time permit several fully analytical manip- 
ulations. We consider the Hernquist (1990) density law: 



P*( r ) = — 



r, + r) 



for which the effective radius R e 
(1983) density law: 



(7) 



1.82 r*, and the Jaffc 



P*(r) = 



M, 



4vr r 2 (r, +rf 



(8) 



for which Re ~ 0.76 r*. 

The density law appropriate for dark haloes in elliptical 
galaxies is yet to be determined. The same kind of density 
profile - though with different masses and scale lengths, i.e., 
with Md and ro replacing M» and r» , respectively - may ap- 
ply to describe the luminous and dark matter distributions 
in a scenario in which both are collisionless and have under- 
gone similar dynamical processes during galaxy formation, 
yet starting from different initial conditions (e.g., Bertin, 
Saglia & Stiavelli 1992, hereafter BSS). On the other hand, 
the distributions of the dark and bright matter may have 
different shapes to the extent that the baryonic component 
has dissipated, thus sinking deeper into the potential well. 

In the first option, DM haloes present a central cusp 
(e.g., Dubinski & Carlberg 1991; BSS; Kochanek 1993, 
1994), and in this mood we investigate HH models, where 
both the luminous and the dark components are described 
by a Hernquist distribution, and JJ models where both fol- 
low a Jaffe profile. In the mood of dissipational collapse, we 
investigate two other classes of models where the dark halo 
density flattens at small radii, while the stellar distribution 
peaks towards the centre. In the first class, that we call HP 
models, the Hernquist luminous component is embedded in 
a Plummer (1911) dark halo: 



Pn(r) 
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(9) 



in the other one (JQ models), the stellar component is de- 
scribed by the Jaffe formula and the dark halo by a trun- 
cated quasi-isothermal distribution: 



po{r) 



' po a rl(rl+r 2 y 



for r < n; 



for r > n. 



(10) 



2.2 The dynamical models 

In order to compare our models with the dynamical prop- 
erties of the observed galaxies, we need their spatial and 
projected velocity dispersion profiles. These are obtained by 
solving the associated Jeans equation (see, e.g., Binney & 
Tremaine 1987, hereafter BT): 



dp*(rK 2 (r) 2a(r)p,(r)a 2 (r) 



dr 



+ 



GM(r) 



with the boundary condition p,(r)a 2 (r) — > for r — > oo, 
where M(r) is the mass within r. We use for a(r) the 
Osipkov-Mcrritt formula: 



a(r) 



1 



a?(r) 



(12) 



(Osipkov 1979; Merritt 1985a,b), so that the velocity dis- 
persion tensor is nearly isotropic inside r a and radially 
anisotropic outside, consistently with N-body simulations 
(see, e.g., van Albada 1982). Analytical expressions for the 
radial velocity dispersion profiles of HH, HP and JJ models 
are given in the Appendix. Their projection is then obtained 
by (e.g., BT, p. 208): 
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dr, (13) 



where £*(i?) = T, I(R) is the surface stellar mass density. 
The solution of (11) then provides the projected velocity 
dispersion profile, via equation (13). 

Having fixed T* over all the FP, a galaxy model is there- 
fore specified by 5 more parameters, namely M*, r» (or R e ), 
Mo, ro, and the anisotropy radius r a . In the following, we re- 
place Mo and ro with the dimensionless ratios 1Z = Mo /M, 
and (3 = ro /r, . 

The observed a entering into the definitions of the fc- 
s in (1) does not correspond to <rp(0), but rather to the 
average over the aperture used for the spectrographic obser- 
vations. That used by D87 for constructing the fc-s of the 
Virgo galaxies in Fig. 1, was normalized to a 4" x 4" aper- 
ture at Coma distance. Being the average of ap(R) over such 
a rectangular aperture very similar to that over a circular 
aperture of 2". 2 radius, we simulate a a by: 



,(-Ra 



2tt 



M*(i? a 



E.(i?) ap(R) RdR, 



(14) 



where Mp(i? ap ) is the projected stellar mass inside 7? ap . We 
therefore mimics the actual observations considering <T ap = 
(j , and we correspondingly get the virial coefficient c% from 
equation (5). The fc-s are then obtained from (3), (4) and 
(6). 

Instead of a fixed angular aperture (as used by D87), 
when calculating cr ap we have for simplicity adopted a fixed 
linear aperture i? ap = 0.02 _R C for all our models. As dis- 
cussed in Section 4.3, such a choice has a negligible effect on 
our results. 



2.3 An analytical approximation 

The described procedure to determine c\ and the fc-s for 
a specified set of parameters r a , TL, and 13 was performed 
numerically solving (11)-(14) for a three-dimensional grid: 
(r a ,ft,/3) € [(r a )min,oo[x[0, 10] x [0.1,40], where (r a ) min is 
the lower acceptable limit for the anisotropy radius as dis- 
cussed in Section 3. We express the resulting c\ in the form: 



A(R, P ) 



C2 - Wm)' (15) 

where A is the virial coefficient in absence of both DM and 
anisotropy, and O represents the correcting factor when such 
ingredients are included. In units of Mqs -2 kpc" 1 km~ 2 the 
value of A(0.02R C ) is 1.74 10 6 and 6.49 10 5 for an Hern- 
P*( r )i (11) quist and a Jaffe stellar distribution, respectively. A good 
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fit (within 5 per cent) for the numerical values of 6 for HH, 
HP and J J models is given by: 

B 



1 + 



-TZ 



C 



1 



E{fi) 



[(rJR e ) + F{0)YV» 



(16) 



with all the coefficients and exponents being reported in 
Tables 1 and 2. 

The adopted form of the function <d retains the main 
physical constraints of the problem. Indeed, — » 1 for r a — > 
co, and TZ —> or (3 — > oo. The linear dependence on TZ 
derives directly from (11). We did not attempt an analytic 
fit for the JQ models. 

2.4 Constraining models to the Fundamental 
Plane 

In order to constrain the models to lie on the FP reproduc- 
ing its tilt, the value of Q at each location on the FP is 
determined using (6), (4) and (15): 

A 



e = 



2ttT» 



10' 



0.26 fci+0.62 



(17) 



and therefore the required trend in either r a , TZ or (3 as 
a function of fci is derived. The stellar mass to light ratio 
Y„ is obtained from (17) with fci — 2.6 and Q — 1, which 
corresponds to assume faintest galaxies to be isotropic and 
devoid of DM. For Hernquist models we find T„ = 5.5, while 
for Jaffe ones we have T» = 2.06. 

3 MAKING THE FP TILT WITH A TREND IN 
THE ANISOTROPY 

In this section we ascribe the entire tilt of the FP to a trend 
with Lb in the anisotropy degree of the galaxies (i.e., in 
r a ), assuming no DM. The values of the anisotropy radius 
at each location on the fci axis are determined by solving 
(17) and then (16) for r a , with TZ = 0. The results are 
shown in Fig. 2 for an Hernquist and a Jaffe stellar dis- 
tribution. The curves are truncated because of the limits 
imposed by dynamical consistency: above a certain lumi- 
nosity, in models constrained to the FP, the phase-space 
distribution function runs into negative values. The lim- 
its can easily be established in the frame of the Osipkov- 
Merritt relation for a(r), without having to know the dis- 
tribution function of the system (Ciotti & Pellegrini 1992). 
We find: r a > 0.25r* ~ 0.1387? e for the Hernquist models, 
and r a >0.05r* ~ 0.036i? o for the Jaffe ones. 

Thus, we conclude that anisotropy alone cannot be at 
the origin of the tilt, because the extreme values of r a that 
would be required correspond to dynamically inconsistent 
models. Note that another argument militate against radial 
anisotropy as the cause of the FP tilt: the requirement of 
radial orbit stability is much more stringent than the simple 
dynamical consistency, and (r a ) m i n increases again. 

4 MAKING THE FP TILT WITH A TREND IN 
EITHER THE DARK MATTER FRACTION OR 
DISTRIBUTION 

Following the negative results of the previous Section, we 
assume global isotropy and move to explore DM as poten- 
tially responsible for the tilt of the FP. We first ascribe all 



Table 1. Numerical values of the fitting /3-indcpendcnt param- 
eters in equation (16). 





B 


b 


C 


c 


D 


d 


HH 


0.076 


1.47 


1.70 


0.96 


0.70 


0.980 


HP 


0.076 


1.47 


1.58 


1.21 


0.95 


1.695 


JJ 


0.009 


1.36 


1.06 


0.48 


0.05 


0.515 



Figure 2. The trend of the anisotropy radius along the FP re- 
quired to produce its tilt, in Hernquist and Jaffe models. The 
curves are truncated at the radius below which the models become 
dynamically inconsistent. The band within dotted lines marks the 
boundaries within which r a can vary at each location on the fci 
axis in accordance with the observed FP tightness. 

the tilt to a trend in the dark to bright mass ratio TZ at con- 
stant f3, and then to a trend in the relative dark and bright 
distributions (3 at constant TZ. 

4.1 Varying the amount of DM 

We set r a = oo, (3 =const, and for 2.6 < fci < 4.4 we 
determine the value of TZ that is required to place the models 
on the FP. Then TZ is obtained from equations (16) and (17) 
for HH, HP and JJ models, and numerically for JQ models. 
Obviously, the larger /3, the larger the variations of TZ that 
are required to produce the tilt. Values of (3<1 may have 
a mere academic interest, although some evidences seem to 
exist in support of a dark halo more centrally concentrated 
than the bright component (Saglia, Bertin & Stiavelli 1992, 
hereafter SBS). By analogy with spiral galaxies, haloes are 
generally considered diffuse (/? > 1), though in some cases 
with significant amounts of DM inside the half-light radius 
(SBS). 

Concerning HH and HP models, for /3 ~ 5 exceedingly 
large values of TZ are required to produce the FP tilt (TZ ~ 
30 — 175), thus wc conclude that an increasing DM content 
from faint to bright galaxies may be at the origin of the 
observed tilt, provided that /3 < 5. 

The same problem affects all the JQ models that we 
have considered, for every values of (3 and r t . 
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Table 2. Some numerical values of the fitting /3-dependent parameters in equation (16) for HH 
(columns 2-4), HP (columns 5-7) and JJ models (columns 8-9; in this case, i(f3) = 12.5 for every 
value of /3 ). 



p 


E(p) 


HP) 


m 


£(/?) 




m 


E(f3) 




0.2 


-1.00 


1.138 


4.40 


-1.00 


1.278 


4.10 


-1.00 


1.143 


0.5 


-1.00 


1.308 


4.70 


-1.00 


1.191 


6.00 


-1.00 


1.235 


1.0 


1.00 


21.67 


4.50 


1.00 


0.946 


3.55 


-1.00 


2.002 


1.5 


1.00 


1.437 


4.60 


2.82 


1.115 


3.73 


1.00 


1.312 


2.0 


1.00 


1.278 


4.40 


3.69 


1.122 


3.58 


1.00 


1.261 


2.5 


1.00 


1.229 


3.90 


3.90 


1.086 


3.39 


1.00 


1.238 


3.0 


1.00 


1.176 


3.75 


4.25 


1.071 


3.28 


1.00 


1.223 


3.5 


1.00 


1.142 


3.55 


4.44 


1.051 


3.16 


1.00 


1.214 


4.0 


1.00 


1.111 


3.45 


4.55 


1.028 


3.07 


1.00 


1.207 


4.5 


1.00 


1.086 


3.35 


4.57 


1.004 


2.97 


1.00 


1.202 


5.0 


1.00 


1.063 


3.30 


4.86 


1.000 


2.92 


1.00 


1.197 


5.5 


1.00 


1.044 


3.25 


5.08 


0.996 


2.88 


1.00 


1.194 


6.0 


1.00 


1.028 


3.20 


4.95 


0.968 


2.79 


1.00 


1.191 



As regards JJ models, 72 never becomes larger than 10 
(for instance 72 ~ 9.5 at the bright end of the FP for (3 = 5), 
thus every value of this parameter is acceptable, for every 
explored value of (3. Fig. 3 (upper panels) show the results 
for HH and J J models, in the cases (3 = 1, 2, 5. 

4.2 Varying the relative concentration of dark 
and bright matter 

We now assume the dark to bright matter ratio 72 to be 
constant among the reference sample of elliptical galaxies, 
and ask the relative concentration of the two components 
(/?) to produce the observed tilt in £3. Thus, the values of 
(3 along the FP are derived for r a — > 00 and 72 — 1, 5, 9. 
For every class of models we find that if 72. = 1, the DM 
in brightest galaxies should be more centrally concentrated 
than the luminous component (/3 < 1). Values of f3 > 1 at 
every location on the FP always require a prevalence of DM 
with respect to bright matter (72 > 1), apart from JQ models 
which are again completely unsatisfactory, their values of (3 
being unrcalistically small for every choice of 72 (we therefore 
reject this class of models). Fig. 3 (lower panels) shows the 
trend of (3 along the FP for HH and JJ models. 

4.3 Aperture effect 

Having used a fixed angular aperture, D87 have sampled 
a larger fraction of the total light (or effective radius) in 
fainter/smaller galaxies compared to brighter/larger galax- 
ies. In fact, the effective radii of the galaxies in Fig. 1 range 
from ~ 0.5 kpc up to ~ 10 kpc, and therefore a circular 
aperture of 2". 2 (i.e., ~ 220 pc radius at the adopted Virgo 
distance), corresponds to a circular region of ~ 0.44 R e ra- 
dius at the faint end of the FP, and of only ~ 0.02 R e radius 
in galaxies at the bright end. When calculating <T ap , we have 
instead adopted a fixed linear aperture of i? ap = 0.02i? c ra- 
dius, thus correctly simulating only the observed a a of the 
brightest galaxies and underestimating the fraction of ef- 
fective radius sampled in the spectroscopic observations of 
fainter galaxies. Thus the derived c|, and the corresponding 
values of the parameter responsible for the FP tilt, are bi- 
ased by such a choice. Indeed, 72 {(3) is set equal to (00) at 
the faint end of the FP, and therefore also the lower limit of 



the driving parameter is not affected by the aperture bias. 
Thus, the derived variation range of the parameters is cor- 
rect, the main effects concerning the intermediate values of 
72 and /3, and the curves in Fig. 3 may be modified in their 
shape only in the range between the starting and the ending 
points. On the other hand, a constant aperture radius for 
all the models have permitted us to express the correcting 
factor 8 in a simple analytical form. 

4.4 Constraints from the tightness of the FP 

The narrow and nearly constant thickness of the galaxies dis- 
tribution about the FP (in the ks direction) corresponds to 
a very small (<12 per cent) dispersion in the ratio of M/Lb 
to the corresponding virial coefficient. If M/Lb ratios and 
virial coefficients are not finely anticorrelated, this implies 
indeed a very small dispersion, separately for both quanti- 
ties, at any location on the FP. In the frame of our basic 
assumption (T* =const), this sets a very severe restriction 
on C2, hence on 6: 

^<0.12, (18) 

which translates into strong constraints on the range that 
each parameter can span at any location on the FP. These 
can be easily derived analytically from equation (16), for the 
three classes of models. For HH and JJ, the dotted lines in 
Fig. 3 represent the band within which galaxy to galaxy vari- 
ations of the corresponding parameter are allowed, and yet 
are consistent with the restrictions imposed by the tightness 
of the FP, i.e., with inequality (18). 

It is evident from these figures that, whatever the struc- 
tural parameter that is responsible for the tilt of the FP, 
and whatever the assumed mass distribution, dramatic fine 
tuning is required to produce the tilt, and yet preserve the 
tightness of the FP (RC). 

Note that also <5r a /r a should be very small at each lo- 
cation on the FP (see the dotted band in Fig. 2), thus once 
more arguing against such an origin of the tilt. 

4.5 Constraints from the velocity dispersion 
profiles 

In the assumption of global isotropy, the observation of 
the radial trend of ap(R) may hopefully give insight on 
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Figure 3. The trend along the FP of the DM content (upper panels) at constant f) and that of the DM concentration (lower panels) at 
constant TL, required to produce the tilt, in HH and JJ models. The band within dotted lines marks the boundaries within which TZ. and 
(i can vary at each location on the k\ axis in accordance with the observed FP tightness. 



the DM content and distribution within the galaxies (e.g., 
SBS; Bertin et al. 1994; Carollo & Danziger 1994a,b; Carollo 
et al. 1995; and references therein). A comparison between 
theoretical and observed op profiles may therefore check the 
reliability of our models and test whether the DM is respon- 
sible for the FP tilt (cfr. RC). In this frame, we have com- 
puted the op profiles of six reference models for every class, 
and for all the explored combinations of TZ and (3. We have 
chosen six FP locations (fcijfe) within the portion of the 
FP actually occupied by Virgo ellipticals (see Fig. 1, open 
circles) : 

Gi = (2.6, 4.2); G 2 = (3.3, 3.2); G 3 = (3.3, 4.2); 
G 4 = (4.2, 3.0); G, = (4.2, 3.6); G 6 = (4.4, 3.4), 

where G = G(ki, fe). G\ corresponds to the faintest model, 
Ge to the brightest one, while models in the pairs G2 , G3 
and G4, G5 only differ for the effective radius (and then sur- 
face brightness). Their luminosities and effective radii are 
determined inverting (1), (2) and (15): 

L B = 2tt i 11Bfc i-°- 62 

^> _ jQl.07fci-0.41 fe 2 -0.21 



and are reported in Tables 3 and 4. For HH, HP and JJ mod- 
els, some representative profiles are shown in Fig. 4 (namely, 
those corresponding to the cases where TZ varies at constant 
(3 = 2, and (3 varies at constant TZ = 5). The values of TZ 
and (3 in each model G; are reported in Tables 3 and 4. 



HH models, especially those at the bright end of the 
FP, are characterized by a sizable central depression in their 
op, while the observed profiles tipically decrease monotoni- 
cally with radius, at least for i?>0.2 R a (Carollo & Danziger 
1994a, b). Therefore the models which better agree with ob- 
servations are those in which the off-centre maximum of op 
lies inside this radius, i.e., those where TZ varies at (3<2, 
and where j3 varies at TZ<5. In the case of HP models, the 
op-profiles always present a prominent off-center maximum, 
and they have to be rejected. On the contrary, the velocity 
dispersion profiles of JJ models are monotonically decreas- 
ing with radius for every explored values of TZ and (3, and 
therefore are consistent with observations. To permit a more 
quantitative comparison between the models and the ob- 
served velocity dispersion profiles, in the JJ case we define 
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Figure 4. Projected velocity dispersion profiles of HH, HP and J J isotropic models with variable TZ and constant (3 = 2 (left panels), 
and with variable P and constant TZ = 5 (right panels). The solid line corresponds to the faintest reference model (Gi), while dotted, 
dashed, long dashed, dot-dashed, and dot-long dashed lines correspond to models G2, Gg, respectively. 



Table 3. HH and HP models: values of Lb an( i Re for the six reference models G;, 
and the corresponding values of TZ for constant fi = 2, and of f3 for constant TZ = 5 
(columns 3-4 for HH models and 5-6 for HP). 

L B R c TZ /3 TZ p 

[f0 10 M Q ] [kpc] (/3 = 2) (71 = 5) (0 = 2) (11 = 5) 

Gi 0.15 0.49 0.0 00 0.0 00 

G 2 0.96 3.41 1.57 3.94 4.75 2.04 

G3 " 1.33 

G 4 10.55 14.91 4.85 2.06 14.66 1.25 

G 5 " 8.48 

G 5 17.95 13.62 5.85 1.84 17.68 1.14 



a slope of ap(R) as: projected velocity dispersion. Table 5 gives the values of Aa 

1 Y)\ for the six G\ models and for R = 2R C , a typical value for 

Aa(R) = 1 y— , (19) the outermost determinations of ap(R) (Bertin et al. 1994). 

ap ^ ' By these values we can recognize characteristic trends in the 

where crp(0) is the maximum (i.e., the central) value of the slope inside each scenario. 
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Table 4. JJ models: values of Lb and Re for the six 
reference models Gi, and the corresponding values of TZ 
for constant (3 = 2, and of /3 for TZ = 5. 



Table 5. Per cent slope of projected velocity disper- 
sion profile as defined in equation (19) for JJ models. 
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In fact, in the assumption that TZ varies from faint to 
bright galaxies at (3 constant less than unity, Aa(2R e ) mono- 
tonically increases along the FP, i.e., the op-profiles system- 
atically become steeper as galaxy luminosity increases. If 
(3 — 1, instead, the slope of the op-profiles is the same for ev- 
ery model, no matter what is the luminosity, while if (3 > 1, 
it systematically decreases along the FP. This can be easily 
understood as the velocity dispersion profiles reflect the po- 
tential well of the systems. Thus, if (3 < 1 an increasing TZ 
along the FP corresponds to an increase of DM content in 
the inner regions of galaxies, while if /3 > 1 external regions 
are involved. In the first case velocity dispersion increases at 
small radii and then crp profiles steepen, while in the second 
one, effects concern external parts of profiles and thus they 
flatten along the FP. If (3 = 1 instead, the potential well 
only deepens, but does not become narrower nor wider and 
velocity dispersion profiles do not vary their slope. 

In the frame of the second scenario (TZ constant and 
(3 decreasing along the FP) velocity dispersion profiles ini- 
tially flatten and then steepen as galaxies luminosity in- 
creases. From previous assumptions Gi is an isotropic and 
DM lacking galaxy, so the required condition fi(G2) < P{G\) 
is equivalent to have added a DM component in the external 
regions of the G2 model. Consequentely the external veloc- 
ity dispersion increases and ap(R) flattens. Moving towards 
the brigth end of the FP, f3 decreases and DM is more and 
more pushed towards central regions and there is a critical 
value /3 C rit when DM starts affecting the central parts of the 
op (R) profiles, rather than their external wings. Thus crp (R) 
steepens for /3 < /3 C rit- 

In conclusion, observations may in principle check 
whether the FP tilt can be ascribed to a trend of TZ at 
13 =const, or it is caused by a variation of (3 at TZ =const, or 
whether a dynamical origin has to be rejected. In the first 
case it is also possible to determine whether DM in galaxies 
is more or less concentrated than the bright component, or 
if they are distributed in the same way. 

5 MAKING THE FP TILT WITH A TREND IN 
THE SURFACE BRIGHTNESS PROFILE 

Systematic deviations of the ellipticals light distributions 
from the standard R 1 ^ 4 profile may also possibly cause the 
FP tilt (Djorgovski 1995; Hjorth & Madsen 1995). In this 
Section we explore such a possibility through a class of mod- 
els in which the log of the surface brightness is proportional 
to R 1/m and m is allowed to vary with galaxies luminosities. 
We assume global isotropy and no DM, thus c 2 depends only 
on the stellar density distribution, and we determine which 
variation of m along the FP is required to generate the tilt. 







Gi 


G2, G3 


G4,, G5 


G 6 


K(/3 = 


0.2) 


57% 


62% 


67% 


68% 


K(13 = 


0.5) 


57% 


59% 


61% 


62% 


TZ((3 = 


1.0) 


57% 


57% 


57% 


57% 


K(P = 


2.0) 


57% 


53% 


50% 


50% 


TZ(p = 


5.0) 


57% 


48% 


43% 


42% 


(i(TZ = 


1.0) 


57% 


53% 


61% 


63% 


p{TZ = 


5.0) 


57% 


46% 


46% 


47% 


(i(TZ = 


9.0) 


57% 


45% 


42% 


42% 



5.1 The models 

The surface brightness distribution of R 1 ^" 1 models is de- 
scribed by the generalized deVaucouleurs law (Sersic 1968): 



I(R) = Jo exp[-b(m)(R/R c ) 1/m ] 



(20) 



where b(m) ~ 2m - 0.324 for 0.5 < m < 10 (Ciotti 1991). 
The dynamical properties of this class of models are deter- 
mined by solving equations (11)-(14) with a(r) = 0. In this 
case the virial coefficient c% depends on the aperture radius 
R ap in a non trivial way. Thus, to correctly simulate real ob- 
servations, equation (14) is solved by averaging crp(R) over a 
fixed angular aperture of l."6 radius, i.e., a suitably varying 
Rap/ 'Rc with galaxy luminosity. 



5.2 The tilt and the tightness 

Assuming faintest galaxies to be R 1 ^ 4 systems, by analogy 
with Section 2, we set: 



C2 



At (flap 

9(m) 



(21) 



where A4 is the virial coefficient for m = 4 and O represents 
the correcting factor when m / 4. In order to produce the 
tilt, this has to increase from 1 to ~ 3 along the FP. The 
values of m that force the six reference models Gi to lie on 
the FP have been correspondingly determined: the result is 
that m has to increase from 4 (Gi) up to ~ 10 (Ge) along 
the fei, as shown in Fig. 5 (upper panel). The figure also 
shows the band within which m can vary for fixed luminos- 
ity consistently with the tightness of the FP. Once again, 
a fine tuning of the driving parameter is required to fit the 
observations: a very small (<10 per cent) scatter of m at any 
location on the FP should be associated to a large variation 
of it with galaxies luminosities. If one assume instead m — 2 
for the faintest galaxies (model Gi), the required variation 
is even larger, about a factor 4, up to ~ 8 for the model 
Gti, and the permitted variation of it at each FP location 
remains very small (Fig. 5, lower panel). 

5.3 Comparison with observations 

It actually turns out that the surface brightness distribution 
of ellipticals is well described by R}^ m profiles with variable 
m, any model with 3 < m < 10 being hardly distinguishable 
from the R 1 ' 4 law in the radial range usually covered by ob- 
servations (Makino, Akiyama & Sugimoto 1990). However, a 
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Figure 5. The values of m for the six reference models Gi re- 
quired to produce the FP tilt in R}/ m models, and the band 
within which m can vary at each location on the fci axis in ac- 
cordance with the observed FP tightness. In the upper panels, 
the faintest model is characterized by m = 4; in the lower panel, 
m = 2 at the faint end of the FP. 

systematic trend of m with galaxy luminosity has recentely 
been reported (Caon, Capaccioli & D'Onofrio 1993, here- 
after CCD), with m increasing from ~ 1 up to ~ 15, thus 
spanning a much wider range than required to produce the 
tilt. Indeed, if one restricts to the Virgo galaxies in com- 
mon with BBF but three (NGC 4406, NGC 4552 and NGC 
4621), m ranges between ~ 2 and ~ 8, in good agreement 
with the required increase of in from faint to brigth galax- 
ies. Worrysome is the apparently large dispersion inferred by 
observations, with m varying by a factor ~ 3 at any given 
luminosity (see Fig. 6), at variance with the observed FP 
tightness. 

At variance with m as possible cause of the FP tilt 
seems a conclusion that may be implicit in the BBF study. 
Being the ratio between tidal radius r t and core radius 
r c about 100-300 for giant ellipticals, when described by 
King (1966) models, one may suppose faintest galaxies have 
r t /r c = 100 and brightest ones r%/r c = 300, thus considering 
a trend in the bright matter distribution along the FP. How- 
ever, as shown in BBF Fig. 5, the corresponding decrease in 
the value of is not sufficient to account for the FP tilt, 
in contrast with our result. However, King models are char- 
acterized by a flat core, at variance with high-resolution 
ground-based and HST observations which suggest an in- 
creasing density towards the very central regions of ellip- 
tical galaxies (Lauer et al. 1992a,b). Therefore we cannot 
consider our result injured by such an argument. 

We conclude that further observational studies are re- 
quired in order to determine whether a progression of light- 



Figure 6. Values of m as fitted by CCD along the major (upper 
panel), equivalent (middle panel) and minor (lower panel) axis 
light profile for the Virgo ellipticals in common with BBF. The 
solid line is the data points best fit line; the dotted lines mark the 
boundary of the permitted variation band of m, in accordance 
with the observed FP tightness (the same as in Fig. 5, lower 
panel). 

profile shapes along the FP really exists among cluster el- 
lipticals. 

6 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated possible structural or 
dynamical origins for the observed tilt of the fundamental 
plane of elliptical galaxies, considering in turn a systematic 
variation along the FP in the radial orbital anisotropy, in 
the dark matter content or distribution, and in the shape 
of the surface brightness profile. In doing so we have varied 
one such parameter at a time, while keeping the other three 
constant. 

Our exploration indicates that all structural/dynamical 
solutions to the fundamental plane problem are rather un- 
appealing, though some are more so than others. This comes 
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from the strong fine tuning that is required, no matter 
whether the driving parameter is the anisotropy radius (r a ), 
the amount of dark matter (71.), its distribution relative to 
the bright matter (/3), or the shape of the surface brightness 
distribution (m). 

In addition to this, we have excluded a trend in the 
anisotropy as possible cause of the tilt because it leads to 
physically inconsistent models, and specific arguments also 
militate against global dark matter content. To produce the 
tilt, the dark to bright matter ratio TZ should increase along 
the FP, from its faint to its bright end. This is just the 
opposite trend that one expects from galactic wind forma- 
tion models (e.g., Arimoto & Yoshii 1987), the only ones so 
far that naturally account for the increase of metallicity (as 
measured by either broad band colors or the Mg2 index) with 
the depth of the potential well (as approximatively measured 
by cr ). Here, the deeper the potential well, the less bary- 
onic material is expelled in a supernova driven wind, thus 
leading to lower final value of TZ. Dissipationless merging 
models - that do not account for the metallicity-<7 correla- 
tion - would predict TZ to remain constant after a merging 
event, or in case decrease slightly, as less bound, preferen- 
tially dark material may escape from the system during the 
merging event. In conclusion, we do not see any good reason 
why the dark matter fraction should systematically increase 
along the FP, and we actually have hints it may decrease 
somewhat. We are therefore inclined to exclude the param- 
eter TZ from being responsible for the tilt. 

Rather more attractive is instead the possibility of a sys- 
tematic decrease of /3 along the FP. Qualitatively, a trend of 
this kind is indeed expected if bright /baryonic matter has 
dissipated deeper into the potential well of small/faint galax- 
ies compared to big/bright ones. Actually, towards the faint 
end of the FP galaxies are characterized by a higher sur- 
face brightness and stellar density, and lower effective radii: 
this suggests that in these galaxies the stellar component 
is more centrally concentrated relative to dark matter com- 
pared to galaxies at the other end of the FP (Guzman, Lucey 
& Bower 1993). Thus, a systematic variation of f3 appears 
more plausible than the previous two alternatives, though 
the fine tuning problem remains. 

Somewhat analogous is the case of a systematic trend 
in the shape of the stellar distribution, for which there ap- 
pears to be some observational support (CCD). If the surface 
brightness profile of ellipticals is well described by a general- 
ized de Vaucouleurs law (7? 1 ^ m -law), then an increase of m 
by a factor of ~ 2 — 4 along the FP is sufficient to produce 
the tilt. The only embarrassment we see with this solution 
is, again, the required fine tuning. 

There remains the possibility of an hybrid origin of the 
tilt, with more than one effect contributing to tilting the FP. 
For example, a small progression of anisotropy, DM concen- 
tration, and shape (m), coupled with a stellar population 
effect causing a modest increase of T*. This is perhaps a 
reasonable solution of the tilt problem, yet a very difficult 
one to test observationally, given that each effect may in- 
dividually be buried in the observational noise, and yet the 
combination of all of them may conjure to produce the ob- 
served tilt. 

Whatever the solution, the tightness of the distribution 
of Virgo and Coma ellipticals about the FP is clear evidence 
for a very standardized and syncronized production of ellip- 



tical galaxies, at least those in clusters. Hypothetical forma- 
tion processes that contain a great deal of stocasticity - such 
as e.g., late merging of spiral galaxies - are likely to generate 
disparate final structures (i.e., r a , TZ, and (5 distributions), 
and stellar age distributions, hence large dispersions about 
the FP. Such scenarios are clearly disfavored by the very 
existence of a tight FP correlation. 

This study has also shown that models where both dark 
and bright components follow a Jaffe density distribution 
(JJ models) may offer a better description of elliptical galax- 
ies. In general, this applies to models where the DM distri- 
bution is similar to the stellar one, albeit less concentrated 
(see also Dubinski & Carlberg 1991; BSS; Kochanek 1993, 
1994). Instead, centrally flat DM distributions frequently 
give physically or astrophysically unacceptable results, such 
as velocity dispersion profiles that steeply increase outward, 
and negative values of their distribution function (Ciotti 

6 Pellegrini 1992). Furthermore, a core radius for the DM 
haloes is hardly justifiable for dissipationless formation, as 
the gravitational force is not characterized by any specific 
scale length. 
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APPENDIX A: VELOCITY DISPERSION PROFILES FOR HH, HP, AND JJ MODELS 



With the assumed radial trend for the of the velocity dispersion tensor, the solution of the Jeans equation can be written in 
integral form, as shown by Binney (1980). Moreover, for our two-component models the solutions can be obtained in closed 
form, and we give here the resulting radial velocity dispersions. The tangential component can be successively obtained by 
using eq. (12). In order to avoid a cumbersome notation, we define the dimensionless variables s = r/r* and s a = r a /r». It is 
straightforward to show that the general velocity dispersion for our models looks like: 

GUj (A„ + s 2 J«.) + TZ(A tU + sj hp) 



p*(r)a T (r) 



(Al) 



27T7-I S 2 + Si 

where the functions A** and 7,, depend on the stellar density distribution only, and 7*d and A,d are the interaction terms 
due to the DM gravitational field. Note that when the velocity dispersion is completely isotropic (i.e. s a — > oo) only the J's 
functions remains; on the contrary, in the formal case of completely radial orbits (s a = 0), the velocity dispersion is described 
by the A's functions. For the Hernquist and Jaffe stellar profiles, the two components due to the pure stellar distribution are 
respectively given by: 

4s + 1 
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We give now the interactions terms for the various distributions. Starting with the HH models, we have: 
1 P+ 1 
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For HP models the functions are more complicated, due to the different qualitative behaviour of the density distributions. 
After lenghty calculations, we find: 

5(2/3 2 - 1) s 3 (10/3 2 - 5) + s 2 (2(3 4 + 9/3 2 - 8) + s(8/3 4 - 9/3 2 - 2) + (5(3* - 10/3 2 ) 



2(1 + /3 2 ) 3 
2/3 4 - lf/3 2 +2 
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Finally, for JJ models, we have: 
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Note that the singularity for /3 = 1 in A*d and 7»d for HH and JJ models is eliminable, for in the limit /3 — » 1 these expressions 
coincide with the corresponding A„„ and /„„. 
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